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Abstract 



In case of superconductors whose electrons attract each other only if they are 

near certain centers, the question arises 'How many such centers are needed to 

make the ground state superconducting?' We shall examine it in the context of 

a random U Hubbard model. In short we study the case where Ui is —\U\ and 

with probability c and 1 — c respectively on a lattice whose sites are labelled 

i using the Gorkov decoupling and the Coherent Potential Approximation 

(CPA). We argue that for this model there is a critical concentration cq below 

which the system is not a superconductor. 
Pacs. 74.20.-z, 74.25.-q, 74.40.-k 
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I. INTRODUCTION 

In many attempts to construct a viable model for High Temperature Superconductors 
the notion of negative - U centers is invokedliTQ. In this connection there is a simple, natural 
question that arises: How many such centers are needed to make a superconductor. In this 
contribution we shall argue that under certain circumstances there is a critical concentration 
Co below which there is no superconducting order. Moreover, we developed a strategy for 
investigating the factors which determine Cq. 

In order to deal with a well posed problem we shall study a Random - U Hubbard Model 
defined by the Hamiltonian 

li — — 2_^ '^ij'^icr'^jcr + 2^ '-^ i'^ia'^icr'^i-a'^i-cr ~ f^ /-^ ^icr^i^y i \^) 

ija i(T icr 

where the coupling constant 

{—\U\ with probality c 
with probability 1 — c 

The question we shall ask is: Is there a finite concentration cq such that for c < Cq the 
cofigurationally averaged, superconducting long range order parameter x vanishes even at 
zero temperature? 

As is natural we define x by the relation 

X=J^Y.X^, (2) 

i 

where the local pairing amplitude is given by 

Xi=< Ci^Cii > . (3) 

Following the conventional notation the bracket < ... > denotes a thermodynamic average 
and 6, for an arbitrary operator O, implies the average of 6 over all configurations Ui, such 
that the fraction of negative U sites is c, with equal weight. A sample will be said to be 
superconducting if x 7^ This implies that x* 7^ on a finite fraction of all sites. Namely, 



if Xi 7^ only on a finite number of sites x will go to zero as A^ ^ oo and the system will 
be regarded as not superconducting. 

To make progress we calculate Xi within the Hartree - Fock - Gorkov (HFG) Decoupling 
scheme for the Greens functions and the averaging over the [/-configurations is accomplished 
with the help of the Coherent Potential Approximation (CPA)0. In short, at the risk of 
missing some important feature of the problem, like localization of electrons, we develop 
a mean field theory for the phenomena described by H in Eq. (|l]) This approach may be 
justified by noting that little is known about the problem at hand systematicallylir a and 
hence as a preliminary study a mean field theory is called for. 

Note that the simplest approximations to the problem would be to set Ui at each site 
equal to its average value U = cU. In some contexts this is called the Virtual Crystal 
ApproximationQ'0. Since, as is well knowncl, any amount of attraction leads to superconduc- 
tivity, U = cU implies superconductivity for all non zero concentrations with the transition 
temperature Tc decreasing albeit non-analytically, with c. Thus before setting out the de- 
tails of the above theory it is worthwhile to pause, briefiy, to consider a number of fairly 
general arguments which suggest that the above conclusion is premature and that there is 
a critical concentration cq of negative U - centers for superconductivity. 



i. Classical Percolation Theory for a mixture of two metals with resistivities pi and p2 
have been studied in the Effective Medium ApproximationETEi. For pi = po and P2 = (Fig. 
1), namely in the case where metal 2 is a superconductor, it yields an effective resistivity 
given by 

1 



Peff 



Po(l — dc) for c < Co 

(Jj 



for C > Co = -; 
a 



(4) 



where d is the number of spatial dimension in which percolation is allowed. Thus, this 



model predicts critical concentrations in c? = 1,2 and 3 dimensions, More over, cq depends 
on the dimensionality d. More generally pe// ~ (c — Cq)** near cq but a mean field theory 
cannot be expected to deal with the critical exponent s adequately. 



ii. The propagation of Cooper pairs between negative —U centers by hopping from 

site to site, where [/ = on the intermediate sites, is depicted in Fig. 2. Assuming that 

the distance between two negative-?/ centers is c~^''^, in units of the lattice constant a on 

a d dimensional lattice, we estimate the number of individual hops / necessary to reach one 

such center from its nearest neighbors. Assuming random walk c~^^'^ = l^/^^^^') (for d > 1). 

If each hop takes h/W seconds where W is the bandwidth for the Cooper pairs, the time 

to travel between two negative U centers is given by r = (h/W) c~^'^~^^^'^. Now we note 

that in between two centers the Cooper pair is without its binding energy U . Consequently, 

such travel is allowed only for such times 5t that the energy uncertainty 5E = h/6t > U . 

Taking (5t = r we conclude that for 5E = Wc^^~^^^'^ > U the pair will propagate for 

SE = Wc^'^^^^^'^ < U the pair will not propagate. Thus for c < Cq, where 

d 



a system of negative U centers will not be superconducting. Presumably, the Cooper pairs 
will be localized. On the other hand for c > cq it will be a superconductor. 



iii. Localization of Cooper pairs by local charge and order — parameter — phase 
fluctuation is the third argument which we wish to recall briefly. It was explored in 
the present context by Doniach and Inula. In a Ginzburg Landau theory on a lattice the 
relative phases of the local order parameters ■j/'j = |'?/'j|e^®' are determined by the quadratic 
term in the free energy function F{{ipi]). This may be written in the form of Josephson 
coupling energies F({0,j}) = ^ X^jj -^i^- cos(9j — Oj), where the precise relationship of the 



coefficients to various parameters of the theory need not concern us here. To describe charge 
fluctuations associated with Cooper pairs arriving and leaving a site a charging energy needs 
to be added to F. Because, the local potential is related to the phase by Josepson voltage 
relation Vi = ^6i this has the form of a kinetic energy term. Finally, to recover a microscopic 
description the local phases are treated as quantum mechanical variables by the Hamiltonian 

^ = i^cEf^^) -^i^,Ecos(e,-e,). (6) 

i \ i / ij 

This is a much studied Hamiltonian in connection with granular superconductivity. In 
particular it was investigated by Gosset and GyorffyM in the Hartree - approximation. In 
short they factorized the wave - function as shown below 

^{{e,}) = ll<P^m (7) 

i 

and found the following self-consistent equation for the individual site wave function 0o(0) 



-l'^c(^y-E,coMe-e) 



4„(e) = £0^0(8) , (8) 



where 



;e^®) = |de0*e^®0o(e)=pe'® 

The amplitude p, determined by solving the above equation numerically is shown in Fig. 
1 (in Ref. |iy),as a function of the ratio Ej/Ec (= Josephson energy/charging energy) For 
Ej/Ec < 0.125 we find p = and hence we conclude that the system of point supercon- 
ductors we have been considering do not have long range superconducting order. Clearly, 
it is tempting to associate Ej with the coupling between the negative - U centers in our 
Hubbard model and assume that it goes to zero as c — *> 0. Evidently, this would imply a 
critical concentration determined by Ec = Ej{cq) . In short, charge fluctuations can de- 
stroy the phase coherence of superconducting order parameter if the coupling between the 
negative - U centers drops below certain critical value. Indeed this was one of the main 
point of the paper by Doniach and Inuicl. In what follows we shall develop a strategy for 



investigating the link between the microscopic model defined by Eq. (Q) and the above 
semi-phenomenological arguments. 

In concluding this introduction we note that the specific task we shall undertake is a 
contribution to the general problem of treating disorder and electron-electron interactions 
simultaneously For a comprehensive discussion of the relevant issues in this field the reader 
is referred to the relatively recent review article by Belitz and KirkpatrickEI. 

II. THE COHERENT POTENTIAL APPROXIMATION FOR THE RANDOM f/ 

HUBBARD MODEL 

The physics described by this simple model appears to be exceedingly rich. For in- 
stance, one might expect that, under some circumstances, the Cooper pairs are subject to 
Anderson localizationO and hence they form a random set of Andreev scatterers for the 
quasi-particlesLJ. Such system of scattering centers may then Anderson localize the quasi- 
particles themselves and turn the system into an insulator below the critical concentration 
Co for superconductivity. However very little systematic fully microscopic work has been 
done on the problem and hence, as a preliminary exercise, a mean field theoretic treatment 
is called for even at the risk of failing to capture some of its important features. In any case, 
as we shall show, even such limited description turns out to be of physical interest. 

Formally, the task is to find the Greens function 



G(z,j;r;{f/,}) 



< T{Qt(r)c+(0)} > < T{QT(r)ca(0)} > 

< T{c+(r)c+(0)} > < T{Q^(r)Qi(0)} > 



(9) 



where the creation and annihilation operators Cjo-(t) and cf^ir) evolve in complex time r 
according the random-f/ Hamiltonian H in Eq. ([l|), T is the r-ordering operator, brackets 
< ... > denote here the usual equilibrium thermal averages corresponding to H, and average 
the result with respect to all arrangement of the [/-centers each denoted by {f/j}. In short 
we wish to find 



{Ui} 

where the probabihty distribution is assumed to be of form 



pi{u^}) = Y[pm 

i 

I c ioT Ui = U 
with P{Ui) = I 

I 1 - c for Ui = 
Note that the local order parameter defined by Eq. (j^) is given by 



(10) 



(11) 
(12) 



(13) 



and hence the knowledge of the averaged one particle Greens functions matrix is sufficient 
to address the question whether or not there is superconducting long range order at a given 
concentration c. 

As we have indicated above we shall now proceed to a mean field approximation to the 
above problem. This consists of two steps. Firstly, we make use of the Hartree-Fock-Gorkov 
decoupling scheme to find the following 'mean-field' equation of motion 



E 



(t^„ + /i - \Uini)5ii + til ^i^ii 
A*Sa (un - /i + ^Uini)5u - ta 



G{l,j]UJn) = IS, 



I] 1 



(14) 



where 



P n 

X^ = ^Ee'""'G'i2(^,^;a;„), 

P n 

\ = -UiXi ■ 

Secondly, we find average of the solution to Eq. (|13]), namely G{i,j; Un', {ui}), over all 
[/-configurations using the Coherent Potential Approximation (CPA). The justification for 
this second step is that the CPA is well known to be a reliable mean-field theory of disorder 
for wave propagation in a medium described by independent random variablesETHll^. 



To implement the CPA we rewrite Eq. (0) in the Dyson form 



G(i,j;cj„) = G°(2,j;u;„) + ^G°(i,/;u;„) Vi G(/,j;u;„) , 



(15) 



where 



V, 



^ \Uini -Ai ^ 



V 



-Af -Wim 



(16) 



The CPA recipe for G{i,j]UJn) is to set it equal to the coherent Greens function 
G'"{i,j;uJn) which is the solution of Eq. ([ISD for the case where the random potential 
Vi is replaced by the energy dependent, complex coherent potential S(c<j„), the same on 
every site. To determine the coherent potential (self- energy) we study, in turn a f/j = —\U\ 
impurity in the coherent lattice. On the impurity site at i we find 



G"(i,i:uJn 



-1 



■iC 



l-G^(i,z;cj„)Vf -E(cu„) G^(z,z;w„), for a = and ?7 , (17) 



where 



/ 1 



ra=0 



and V 



a=U 



Um 



A, 



V 



-A* - |f/n, 



(18) 



Then, the usual CPA condition which determines the self-energy S(u;„,) is given by 



:G(°) 



Ctx^ '[t,i;UJn_ 



+ il-c)G^^\t,i;iUn) = G^it,i;iUr, 



(19) 



Similar equations have been used to describe random superconductors by LustfeldEJ and 
more recently by Litak et a£3. The principle difference between our present concerns and 
that of these earlier authors is that we are focusing on the randomness of the interaction 
parameter Ui and not on the random site energies Cj as was their aim. To put it on other 
way we are studying a problem analogues to that of ' spin-glass '□ rather than that of dirty 
superconductors. 

Equations (PJH) and Eq. (W^) together with 



P n 

x" = |Ee'""'G?2 (20) 

P n 

n = cn^^^ + (1 - c)n(°) 
Y = cx^u) + (I - c)x^^\ 

where a = and U as before (Eq. [1^), are the fundamental equations of our theory. 
Manipulating the CPA equations yield the following gap equation: 



P n 



G^^^W ,.„nn.. ( 2cj„-TrS||G(^)| 



.Tr{G(^)}i^ + ||G(^) 



c- 



^^n ' \\G \\ \ ^^n \\G 



X ■ 
(21) 



In what follows we present results of solving the above equations numerically for various 
interesting regimes. Of particular interest is the large U limit. As Ui change its values form 
to — |f/| there exists Mott-Hubbard metal-insulator transition for large enough interaction 
\U\. An other interesting feature of the problem at hand is that fluctuations of pairing 
potentials Aj, which changes randomly from to A^, invalidate the Anderson theoremt-3ll3 
and hence states appear in the gapE2l. 

III. ORDER PARAMETER FLUCTUATIONS 

At first we have calculated Tc and x for ^sro temperature (T = 0) by means of VGA 
where, as we mentioned in the introduction effective interaction between electrons f/e// = 
cU . Figs. 3a and 36 show the critical temperature Tq normalized to the corresponding 
quantities of the clean system with U on every site, namely Tc{c= 1) and averaged pairing 
parameter A(T = 0, c = 1), calculated for effective interaction U^ff respectively as functions 

9 



of concentration c. Calculations were done for various values of the interaction parameter 
U/W: -0.3, -0.5, -1.0, -2.0 and for a half filled band: n = 1. One can see that for these 
approximations there is no evidence of critical concentration cq > below which the systems 
is normal at T = i.e. no percolation. 

As is clear from Eq. (0), f/j fluctuating between and ~\U\ has two distinct direct 
consequences. On the one hand it causes the Hartree potential ^UiUi to fluctuate. On 
the other it gives rise to a fluctuating pairing parameter - Aj. As it turns out these two 
effects have very different influence on the solutions to Eqs. ( |T8i|l9tp0i^Il) . Therefore, we 



examine them separately. As disorder was treated by CPA, at first we made calculations 
after neglecting Hartee potential ^UiUi, in Eq. (|18D and studied the case of order parameter 



fluctuation on their own. This means that we took the impurity potential in Eq. (|T8D to be 

( -A,^ 
V, = . (22) 

1^-Ar ) 

In Figures 4a, b, c we show the critical temperature Tc/Tc{c = 1) (a), the order param- 
eter (for T = 0) x/x(c = 1) (b) and the local pairing potential on U site Aij/Au{c = 1) (c), 
versus concentration of negative centers c for n = 1. Calculations were done by means of 
CPA neglecting the Hartee term and using the same values of interaction parameters as in 
Fig. 3. 

Surprisingly, our simplified CPA results agree with the VCA argument in the introduc- 
tion in as much as we found non-zero local order parameter on both the U = and U < 
sites at all concentrations c 7^ 0. That is to say we obtained finite x 7^ and Tc for any 
value of concentration c and interaction f/j < and no evidence of percolation. The order 
parameter x increases gradually from to its maximal value with changing the concentra- 
tion c. Interestingly, in the large U limit, |f/|/14^ > 0.5, Tc and Ajj are nearly constant for 
various c and they reach large finite values for arbitrary small concentrations of negative 
centers c. 

To return to the full CPA solution we have used the full impurity potential \/"=^ as in 
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Eq. ([T8|) . Figures 5a, b, c, d show the critical temperature T^' /Tc{c = 1) (a), the order 
parameter x/x(c = 1) (b) and the local pairing potential on U site Au/Aij{c = 1) (c), 
versus concentration of negative centers c for n = 1 and the same interactions as in Figs 3 
and 4. Here one can clearly see that all that quantities Tc, x{T = 0) ^^'^ '^u are tending 
to zero for some small enough concentration of negative U centers Cq. Below this critical 
concentration the system is normal. For larger interaction {U/W = -1.0, -2.0) cq = 0.5 and 
the order parameter scales as x ~ (c — Cq)^/^. Decreasing \U\ we observe systematic decrease 
of Co. 

To investigate the case of critical concentration cq we have studied the density of quasi- 
particle states both in superconducting and in the normal states. In the latter case, for 
large enough interaction \U\ > 0.5W, there exists a band splitting in the system. With 
changing the concentration c we observe Mott metal-insulator transition. It is caused by 
large fluctuations of Hartree term ^UiUi (Eq. [M]) as in the original paper of HubbardEj. In 
Figs. 6 a, b, c we plotted the densities of states (full line) and the local density of states 
on U site (dashed line) for U/W = —2.0, n = 1 and T = for different concentrations c: 
c = 0.4 (a - normal metal), c = 0.5 (6 - insulator) and c = 0.6 (c - superconductor). The 
Fermi energy in these plots: ep = fi = 0. Thus changing c from to 1 system changes from 
normal metal (Fig. 6a) to a superconductor (Fig. 6c) through an insulator (Fig. 66). Re- 
markably, for a low concentration of negative centers (Fig. 6a) c = 0.4 < cq (here Cq = 0.5) 
in spite of finite and relatively large value of averaged density of states at the Fermi energy: 
D(0) = -ilmGfi(0 + i5), the local density of states on U sites D^{0) = -^^^^^(O + i5) 
appears to be extremely small. Evidently the doubly occupied states form a lower 'Hubbard' 
band split off from the upper band which is associated with the singly occupied sites. 

This effect has been further investigated for other band fillings. The transition from 
normal to superconducting phase occurred for each of band fillings n at some, specific, 
critical concentration Co(n). Figures 7a, b, c, d show simultaneously the order parameter 
X (a), the local pairing potential on U site Au (b), the local charge on U site riu (c) and 
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the chemical potential /i (d) plotted versus concentration c for U/W = —2.0, at T = 
and several values oi n {n = 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2) Interestingly, that 
transition from superconducting to normal phase is accompanied by a large value of local 
charge occupation riu (Fig. 7c) and large jump of a chemical potential ^ (from one subband 
to another) near cq (Fig. 7d). It appears that for c below cq njj ^ 2. Namely, every U site 
is doubly occupied with a pair of electrons (Fig. 7c). Because there are no empty spare U 
sites in the system these pairs cannot move. That is to say, they are localized on the U sites. 

Similar calculations have been performed for smaller interaction \U\ {U/W = —0.5) The 
corresponding results are presented in Figs. 8 a-d respectively. Here the interaction \U\ 
is not large enough to create a band splitting effect but the tendency with x — > is still 
observable as concentration c is tending to some finite Cq > (c — i> Cq). Here Cq is less that 
in former case of larger interaction \U\ = 2W (Fig. 7). The occupation of negative centers 
is larger than n but clearly less than 2 electrons per site. For small enough band filling 
n = 0.2, the order parameter x was finite for all c and we have not observed a percolation 
phenomenon. For larger band fillings we have obtained x = below c < cq but instead of a 
square root behaviour x ^ {c — cq)^/^ for c close to cq for larger U (U/W=-2.0 fig. 7a) here 
X goes to zero rather in the asymptotic way (Fig 8a). 

To investigate the demise of the superconducting state near cq we have studied the 
density of states in the appropriate region of parameter space. Figures 9a and b shows 
the quasiparticle densities of states (a) and the local densities of states on U site (6) for 
U/W = —0.5, n = 0.4 and three values of c specified in the figures. It is clearly visible 
how the superconducting gap is filled, due to pair breaking, in with c. Beginning from 
the clean system with interaction U in every site c = 1 we start from the sharp edges in 
the quasiparticle density of states (Fig. 9a) then for smaller values of c (c = 0.6) the gap 
parameter A is of the same order (Fig 76) but the real gap in the quasiparticle density of 
states D{E) changes significantly. The gap becomes smaller with smaller c and looses its 
clear edges. For small enough c {c = 0.14) it nearly disappears. Clearly, the Anderson 
theorem for a superconuctor with nonmagnetic disorder is not satisfied in this casetZl. As is 
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well known, according to Anderson theorem the gap remains absolute in presence of disorder 
due to potential scattering provided the spatial fluctuations of Aj about A are negligiblell3. 
Clearly, in the random interaction case this is not true and this kind of disorder leads to 
pair breaking. 

Thus on account for the large fluctuations of pairing potential Aj in our system,due to 
disorder, we observe a qualitative change in quasiparticle density of states shown in Fig 9a. 
These fluctuations lead also to complicated gap equations where Tc is determined not only 
by G^ but also by S, G^ (Eq. P). 

Finally, we investigated the factors which determine the critical concentration cq. In Fig. 
10 we show Co as a function of band filling for two interaction parameters U/W = —2.0 and 
—0.5. In both cases function can be approximated by a straight line cq = a + bn. In case of 
U/W = -2.0, a = and 6 = 0.5 but for U/W = -0.5, a ^ 0.32 and b ^ 0.6. 

IV. CONCLUSIONS 

We have examined the question of percolating superconductivity in the context of a ran- 
dom U Hubbard model. We have studied the case where f/j is —\U\ and with probability 
c and 1 — c respectively on a lattice whose sites are labelled i using the Gorkov decoupling. 
Changing concentration c we checked that simple averaging procedures like Virtual Crystal 
Approximations (VCA) do not lead to any zero temperature phase transition. Furthermore, 
we found that if charge fluctuations are neglected even a full mean fleld theory of disorder, 
like the CPA, does not predict a percolation transition. However, when the fluctuations in 
Hartree potential are included on equal footing with the fluctuations in the pairing poten- 
tial A and the problem is treated in the Coherent Potential Approximation a percolation 
phenomena, with a critical concentration cq of the negative U centres, is discovered in our 
fully microscopic theory. For c < cq the lack of superconductivity is due to Mott localiza- 
tion of Cooper pairs and its high lights the qualitative difference between disorder in the 
crystal potential and the disorder in the interaction between the carriers. Having found the 
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critical concentration cq we investigated its dependence on various parameters which defined 
the problem. In short we studied Co(n, f/). For strong attractive interaction cq = n/2 and 
X ^ {c — Co)^/^ near cq but for smaller interaction x ^ (as c — > cq) rather in a non 
polynomial manner. Calculations have been performed by a real space recursion algorithm 
which we developed for disordered superconductors in earlier publicationE^. 
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FIGURES 

FIG. 1. A mixture of normal, with resistance p = po, and super, with p = 0, conductors. 

FIG. 2. The propagation of Cooper pairs between negative U < centers by hopping from 
site to site, where U = on the intermediate sites. 

FIG. 3. The critical temperature Tc/Tc{c = 1) (a) and the pairing parameter A/A(c = 1) 
(6) (calculated with VGA) versus concentration of negative centers c for n = 1. Values of the 
interaction parameter U specified in the figures. 

FIG. 4. The critical temperature Tc/Tc{c = 1) (a), the order parameter x/x(c = 1) (b) and 
the local pairing potential on U site Au/Au{c = 1) (c), calculated with GPA neglecting diagonal 
Hartree terms, versus concentration of negative centers c for n = 1. Values of the interaction 
parameter U specified in the figures. 

FIG. 5. The critical temperature Tc/Tc{c = 1) (a), the order parameter x/x(c = 1) (b) and 
the local pairing potential on U site Au/Au(c = 1) (c), calculated with CPA including diagonal 
Hartree terms, versus concentration of negative centers c for n = 1. Values of the interaction 
parameters U specified in the figures. 

FIG. 6. The densities of states (full line) and the local density of states on U site (dashed line) 
for U/W = —2.0 and n = 1 for different concentrations c: c = 0.4 (a - metal), c = 0.5 (b - insulator) 
and c = 0.6 (c - superconductor). 

FIG. 7. The order parameter x (a), the local pairing potential on U site Ajj (b) and the 
local charge on U site njj (c) and the chemical potential /i (d) plotted versus concentration c for 
U/W = -2.0 and several values of n (n = 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2 - the direction of 
n changing is pointed out by the arrow). 
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FIG. 8. The order parameter x [o)-, the local pairing potential on U site A(7 (6) and the 
local charge on U site riu (c) and the chemical potential ^ (d) plotted versus concentration c for 
U/W = -0.5 and several values of n (n = 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2 - the direction of 
n changing is pointed out by the arrow). 

FIG. 9. The quasiparticle densities of states (a) and the local densities of states on U site (b) 
for U/W = —0.5, n = 0.4 and several values of c specified in the figures. 

FIG. 10. The critical concentration of negative centers c versus band filling n for U/W = —2.0 
and —0.5. 
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FIG. 1. A mixture of normal, with resistance p = po, and super, with p = 0, conductors. 
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FIG. 2. The propagation of Cooper pairs between negative U < centers by hopping from site 
to site, where U = on the intermediate sites. 
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FIG. 3. The critical temperature TcjTcic = 1) (a) and the pairing parameter A/A(c = 1) 
(6) (calculated with VGA) versus concentration of negative centers c for n = 1. Values of the 
interaction parameter U specified in the figures. 
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FIG. 4. The critical temperature Tc/Tc{c = 1) (a), the order parameter x/x(c = 1) (6) and 
the local pairing potential on U site A[//A[/(c = 1) (c), calculated with CPA neglecting diagonal 
Hartree terms, versus concentration of negative centers c for n = 1. Values of the interaction 
parameter U specified in the figures. 
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FIG. 5. The critical temperature Tc/Tc{c = 1) (a), the order parameter x/x(c = 1) (6) and 
the local pairing potential on U site Au/Au{c = 1) (c), calculated with CPA including diagonal 
Hartree terms, versus concentration of negative centers c for n = 1. Values of the interaction 
parameters U specified in the figures. 




FIG. 6. The densities of states (full line) and the local density of states on U site (dashed line) 
for U/W = —2.0 and n = 1 for different concentrations c: c = 0.4 (a - metal), c = 0.5 (6 - insulator) 
and c = 0.6 (c - superconductor). 
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FIG. 7. The order parameter x i^)-: the local pairing potential on U site A[/ (6) and the 
local charge on U site njj (c) and the chemical potential n (d) plotted versus concentration c for 
U/W = -2.0 and several values of n (n = 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2 - the direction of 
n changing is pointed out by the arrow). 
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FIG. 8. The order parameter x i^)-: the local pairing potential on U site A[/ (6) and the 
local charge on U site njj (c) and the chemical potential n (d) plotted versus concentration c for 
U/W = -0.5 and several values of n (n = 1.8, 1.6, 1.4, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2 - the direction of 
n changing is pointed out by the arrow). 
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FIG. 9. The quasiparticle densities of states (a) and the local densities of states on U site (6) 
for U/W = —0.5, n = 0.4 and several values of c specified in the figures. 
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FIG. 10. The critical concentration of negative centers c versus band filling n for U/W = —2.0 
and —0.5. 



